Mott transitions in ternary flavor mixtures of ultracold fermions on optical lattices 
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Ternary flavor mixtures of ultracold fermionic atoms in an optical lattice are studied in the case 
of equal, repulsive on-site interactions (7 > 0. The corresponding SU(3) invariant Hubbard model is 
solved numerically exactly within dynamical mean-field theory using multigrid Hirsch-Fye quantum 
Monte Carlo simulations. We establish Mott transitions close to integer filling at low temperatures 
and show that the associated signatures in the compressibility and pair occupancy persist to high 
temperatures, i.e., should be accessible to experiments. In addition, we present spectral functions 
and discuss the properties of a novel "semi-compressible" state observed for large U near half filling. 
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Starting with the achievement of Bose Einstein con- 
densation in 1995, the creation and study of quantum 
degenerate atomic gases has led to discoveries with enor- 
mous impact far beyond atomic physics [l[. In particu- 
lar, atomic gases can be driven to the strongly correlated 
regime by switching on optical lattices and/or using Fes- 
hbach resonances. In the bosonic case, the localized Mott 
phase is accessible already in single-component (lattice) 
systems [5] . With the recent observation [sl, ^ of Mott 
transitions in balanced 2-flavor (i.e. 2 hyperfine state) 
mixtures of fermionic '^'^K atoms with repulsive interac- 
tions, such systems are now established as highly tunable 
quantum simulators of condensed matter ^5i]. 

At the same time, ultracold atoms offer new degrees of 
freedom: In contrast to the electronic case, with only two 
spin states, fermionic atoms have large hyperfine multi- 
plets. Early experiments ^6*1 involving three hyperfine 
states of ^^K have prompted theoretical investigations of 
3-flavor mixtures on optical lattices. These studies have 
focussed on the case of pair-wise equal attractive interac- 
tions which may induce trionic phases and exotic super- 
fluidity with a 3-component order parameter, somewhat 
analogous to QCD 0]. Within the last year, balanced 3- 
flavor ^Li mixtures have been trapped and studied across 
Feshbach resonances [1, Soon, such systems, both 
with repulsive and attractive interactions, should be re- 
alized also in optical lattices. 

In this paper, we explore the properties of fermionic 
3-flavor mixtures on optical lattices in the case of repul- 
sive interactions. While the many ordering patterns con- 
ceivable in such systems will certainly warrant extensive 
studies at some point, ordering phenomena have escaped 
experimental detection even in the simpler 2-flavor case 
so far. We will, thus, concentrate on the Mott physics 
of homogeneous phases. Note that this regime is partic- 
ularly challenging since it requires nonperturbative the- 
oretical approaches. 

Ternary mixtures of fermions on an optical lattice can 
be modeled via the Hubbard-type Hamiltonian 

^ — ^ ^ ^ia^ja^ ^ ^ Uaa'^ia^ia' ^ ^ t^a ^ia • (1) 

{ij) ,a i,a<a' i,a 

Here, {ij) denotes nearest-neighbor sites; a G {1,2,3} 



labels the fermionic flavors, t parameterizes the hopping, 
U the on-site Coulomb interaction, and fj, is the chemi- 
cal potential. In Eq. we have neglected higher Bloch 
bands and the confining potential. For fermions of a sin- 
gle species (in the vibrational ground state), ta ~ t^' . 
Following the literature, we will also assume pairwise 
equal interactions and study the SU(3) symmetric limit 
ta = t, Uaa' = U, fia = fJ-+U ; by this definition, particle- 
hole symmetry corresponds to a sign change in fi. 

This system is treated within dynamical mean-field 
theory (DM_FT) which retains the dynamics of local cor- 
relations [10]. The semi-elliptic "Bethe" density of states 
with bandwidth W = At* mimics [ll| a simple cubic lat- 
tice (with t — t*/y/6); as usual in DMFT studies, the 
scale is set by the scaled hopping amplitude t* — \fZt 
for coordination number Z . The DMFT impurity prob- 
lem is solved using the Hirsch-Fye quantum Monte Carlo 
algorithm (HF-QMC) [l^, first in the conventional form, 
with imaginary-time discretization Ar < 0.8/[/; the va- 
lidity of these results is later verified and unbiased results 
are obtained using the numerically exact multigrid imple- 
mentation [l^. Note that already conventional HF-QMC 
is competitive with continuous-time QMC 'T3| ; the multi- 
grid variant is even more efficient. All QMC based meth- 
ods share the advantage, compared to numerical renor- 
malization group (NRG) approaches, of being reliable at 
the experimentally relevant (elevated) temperatures. 

Results at low temperature - Theoretical studies of 
electronic 1-band models are often restricted to half fill- 
ing (/i = 0). Mott metal- insulator transition are, then, 
signaled by kinks or jumps in thermodynamic properties 
(as a function of interaction U and/or temperature T) 
or, more directly, by the opening of a gap in the local 
spectral function A{uj). In the present 3-fiavor case, gen- 
uine Mott physics can be expected only far away from 
particle-hole symmetry; consequently, the chemical po- 
tential fi is an essential additional parameter. In the 
following, we will choose a relatively low temperature 
T = t* /20 and explore the ([/, /i) space with particu- 
lar emphasis on properties that are most accessible in 
quantum gas experiments and which are related to Mott 
physics. Temperature effects and spectral functions are 
to be discussed later. 




FIG. 1: (Color online) HF-QMC estimates of particle density 
n{fj,) atT — t* /20 for various on-site interactions U. Plateaus 
at integer filling indicate localized Mott phases (for U/t* > 6). 
In the inset, curves are shifted (in steps of 0.1) for clarity. 



FIG. 2: (Color online) HF-QMC estimates of pair occupancy 
D at r = t* /20 vs. chemical potential jj, (main panel) and 
vs. density n (inset), respectively. The suppression of D with 
increasing U is clearly enhanced on the Mott plateau at n « 1. 



Figure [T] shows the filling n — "^^{ha) versus the 
chemical potential /i for a range of on-site interactions U . 
Initially, for U — Q, n varies smoothly and rapidly with 
/i from an empty band (n = 0) at /i/i* < —2 to a full 
band {n ~ 3) at ji/t* > 2. With increasing U, the slope 
generally decreases, but the curves remain smooth until, 
for U/t* > 6, plateaus develop at integer filling n — 1, 
n = 2 which signal the onset of localized Mott phases 
(and correspond to gaps in the spectral function). The 
particle- hole symmetry around /i = is evident; from 
now on, we will only consider /i < 0. The plateaus ap- 
pear perfectly fiat even at the magnified scale of the inset 
of Fig. [T] (for n ~ 1); tendencies towards hysteresis are 
seen for U/t* =6 (at n/t* « -2.5). We wiU argue below 
that T = t* /20 is very close to the critical temperature 
T* of a second-order Mott critical end point (cf . Fig. [3]) . 

From our results, detailed predictions for experimen- 
tal measurements can be derived within the local density 
approximation (LDA). For a specific trap potential and 
atom number, the data of Fig. [H then, translates to oc- 
cupancies Ui of the inhomogeneous system, integrated 
column densities and cloud sizes. Analogous approaches 
have been used in the Bloch-Rosch collaboration for es- 
tablishing the Mott phase in the 2-flavor case [1, • 

A fundamental property of Mott phases is the suppres- 
sion of density fiuctuations; in the 2-flavor case these are 
usually quantified by the double occupancy {hi-^hii). In 
the present 3-flavor context, we consider the more gen- 
eral pair occupancy D = X]q<q' ("-M^-ia')' which retains 
the energy relation E = E'kin + UD [16]. As seen in the 
main panel of Fig. [21 D depends strongly on fi and U; 
the dependence on /i is mostly monotonic, except for the 
vicinity of plateaus (for U/t* > 6). The impact of U 
is best understood at fixed density n (see inset of Fig. 
[5]): starting from the noninteracting limit [7 = 0, where 
D = D is suppressed with increasing U at all 

n. This suppression is strongly enhanced at n = 1 for 



U/t* > 6. Fortunately, the probabilities for n-fold oc- 
cupancy can be easily measured experimentally. In fact, 
this technique has been used by the Esslinger group for 
detecting Mott phases in the 2-fiavor case via changes in 
the slopes dD/dn close to integer filling [S]; analogous 
slopes are marked by thick line segments in the inset of 
Fig. [2 

One interesting feature of Fig. [1] not discussed so far 
is that the slopes dn/dfj, seem to saturate in the range 
—2 < fi/t* < 2 for large U. The corresponding com- 
pressibility K = dn/dfjL [17|, computed via numerical dif- 
ferentiation of the data of Fig. [U is shown in Fig. [3l 
Note that - in contrast to expressions such as n^^dn/ dfi 
- this definition of k, retains the particle-hole symmetry 
(k(— /i) = i^ifJ-)) and fulfils a sum rule n{iM)d^ — 3), 

just as the spectral function. Indeed, reduces to 

the slightly broadened semi-elliptic spectral function for 
U — Q (thick solid line in Fig. [3]). With increasing U , 
the height of this central peak is rapidly reduced and 
the weight of K(/i) is shifted towards larger with a 
subpeak at |^| « [/ + t* developing for U/t* > 2. The 
intermediate minima are initially smooth; their positions 
coincide with integer filling as is apparent from the curves 
K{n) shown in the inset of Fig. [3] Finally, for U/t* > 6, 
sharp gaps appear; these Mott plateaus [with K(/i) w 0] 
collapse to points in the representation Hi{n). 

Coming back to the saturation of slopes dn/dfi, let 
us first note that the peak value ^(/i = 0) reaches a fi- 
nite limit of about 0.2/t* for large U (and low T), i.e., 
scales with the inverse of the (only remaining) energy 
scale W oc t* . Even more importantly, the variation of 
K(/i) within the central compressible phase decreases with 
increasing U until k is nearly constant for U/t* > 8 [in 
contrast with the edge phase (n < 1) where Kt* varies 
between and about 0.3]. This additional plateau marks 
a new "semi-compressible" state sandwiched between 2 
equivalent Mott phases which is specific to SU(2Af -f- 1) 
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FIG. 3: (Color online) HF-QMC estimates of compressibility 
K = dn/dn at T = t*/20 as a function of chemical potential 
fj, (main panel) and density n (inset), respectively. 

systems (M > 1) [l^ and should be easy to probe exper- 
imentally. 

Impact of temperature - While experiments with 
fermionic atomic gases in optical lattices offer much 
greater flexibility in varying interaction and hopping pa- 
rameters (compared to solid state experiments), they 
are currently restricted by initial temperatures not much 
lower than the Fermi temperature Tp. On the other 
hand, true Mott transitions and/or ordering phenomena 
typically occur on energy scales that are 1-2 orders of 
magnitude below the kinetic energy scales (such as Tp or 
the bandwidth W). Therefore, it is essential to study the 
systems of interest also at elevated temperatures and, if 
possible, to identify signatures of Mott transitions which 
persist to temperatures much higher than the critical 
temperatures, i.e., signal well-defined crossover regions. 
At the same time, good estimates of the critical temper- 
ature(s) T* of the Mott transitions are desirable. 

In Fig. 21 we present multigrid HF-QMC estimates 
(thick lines) of n{fi) (top panel) and K{fi) (lower panel) 
for a wide range of temperatures at a relatively strong in- 
teraction U/t* = 6. Clearly, the temperature effects are 
small across wide regions of /x, with deviations in n typi- 
cally of the order 10~^ for T < t* /5. However, extended 
hysteresis regions at the plateau edges are observed only 
at the lowest temperature, T = r/30. At T = t* /20, co- 
existence is hardly significant (around fi/t* « —2.5) while 
the compressibility is strongly enhanced in both transi- 
tion regions, as expected very close to 2"^^ order critical 
end points |17| . Consequently, we estimate the critical 
temperature as T* = (0.05 ± 0.01)r. At higher tem- 
peratures, all curves are smooth. Nevertheless, a Mott 
plateau can be identified up to T = t*/5: in the central 
region /i/t* « —3.3, the compressibility is still reduced 
by two orders of magnitudes compared to off-plateau re- 
gions. This is good news for experimentalists: a factor 
of 4 above the critical temperature [roughly correspond- 
ing to an entropy s = log(3) per particle], the (remnant 
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FIG. 4: (Color online) Multigrid HF-QMC estimates of den- 
sity n{fj,) (top panel) and compressibility «:(/^) — dn{fi)/dfi 
(bottom panel) at interaction U /t* = 6 for various tempera- 
tures T. The strong suppression of k at n ~ 1 survives up 
to T ~ t* /5 3> r* ~ t* /20. For comparison, thin solid lines 
indicate HF-QMC data for T = t* /20 (cf. Figs. [T] and O. 

of the) Mott plateau should still be well visible. Even 
at T = t* /2, a full order of magnitude above T*, the 
compressibility is stronger reduced, in the Mott crossover 
region, than at the subcritical interaction U jt* = 5 for 
low T (cf . Fig. [3]) . A similar picture emerges for the pair 
occupancy (not shown): both the plateaus and the non- 
monotonic behavior in D(/x) persist up to T < t*/5. 

Also shown in Fig. [J] are conventional HF-QMC esti- 
mates for T = t*/20 (and Art* = 0.12; thin solid lines). 
The good agreement with the numerically exact multi- 
grid data (thick solid lines) shows that the discretization 
effects of the former method are hardly significant; thus, 
the data shown in Fig. [1]- [3] are reliable. 

Spectra - Finally, let us discuss an experimentally 
more challenging observable: the local spectral func- 
tion A(lu). Figure [5] shows estimates obtained via maxi- 
mum entropy analytic continuation from numerically ex- 
act multigrid HF-QMC Green functions for U/t* = 6. 
At half filling (top curve), A(lu) is symmetric and shows 
a rough 3-peak structure (with peaks at a; = and 
Lu/t* sa ±7), reminiscent of the familiar l-band-2-spin 
correlated metallic phase. However, important differ- 
ences arise in the 3-flavor case: (i) the Hubbard bands, 
although well separated, carry little weight (of about 10% 
each); (ii) the central peak (at \uj/t*\ < 4) has a three- 
fold substructure; (iii) the quasi-particle peak remains 
for arbitrary U (with only 10% reduction in weight for 
5 < U/t* < 10; not shown). 
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FIG. 5: (Color online) Multigrid HF-QMC estimates of 
spectral function A{uj) for T = t* /lO, U/t* = 6 at 
different values of /i, corresponding to densities n = 
1.5, 1.27, 1.07, 1.00,0.9996, and 0.91 (from top to bottom). 

With decreasing filling, triply occupied sites become 
even less important: the upper Hubbard band vanishes 
for fi/t* < —2. In contrast, the lower band rapidly 
acquires weight: about 26% for pi/t* = —1, while the 
quasiparticle peak shrinks at nearly constant peak height 
{A{0)t* « 0.3). It disappears when /i reaches the Mott 
plateau: at fi/t* = —3.2, a well-defined central gap has 
opened. The remaining peaks have centers of masses at 
uj = —2.8t* and oj = 3.2t*, respectively; the difference ex- 
actly equals the interaction U, just as one would expect 
in the half-filled 2-spin case. However, in the present 3- 
fiavor case (also at n = 1), the weights are strongly asym- 



metric: 1/3 and 2/3 for the lower and upper Hubbard 
band (of which the latter has formed from the central 
peak at 71 = 1.5), respectively. For further decreased fi, 
a quasiparticle peak reappears from the lower subband. 

Conclusion - In this first nonpcrturbative study of 
the repulsive SU(3) symmetric Hubbard model, we have 
thoroughly explored its homogeneous phases at low and 
intermediate temperatures in the (U, fi) space. We have 
shown that signatures of Mott transitions close to integer 
filling persist to elevated temperatures, facilitating ex- 
perimental observation of strong correlation phenomena, 
e.g., in balanced 3-flavor mixtures of ^"K (with roughly 
equal positive scattering lengths at small magnetic field 
6]) or ^Li d, 3 on optical lattices; systems with Raman 
excited states 19] provide an interesting alternative. 

We have found peculiar features in spectra and a novel 
"semi-compressible" phase with asymptotically universal 
K (independent of U, T, and fj, or n), which is clearly 
linked with the presence of adjacent equivalent Mott 
lobes. Experimental studies of this new physics might be 
challenging, but are within reach of existing techniques. 

Our work could be extended towards full quantitative 
correspondence with upcoming experiments by explicit 
inclusion of the trap potential, within LDA (similar to 
recent studies [1, [13, [2l[ for 2-flavor systems) or beyond 
; also, ordered phases [23,[13| could be considered. 
It seems, however, more urgent to investigate which prop- 
erties of SU(3) symmetric systems survive for unequal 
interactions (with less severe 3-body losses [i,[i|), when 
flavor-selective Mott physics [l^l may be expected. 
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